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Abstract 

Very high energy (VHE, E >100GeV) y-rays are absorbed via interaction with low-energy photons from the extragalactic back- 
ground light (EBL) if the involved photon energies are above the threshold for electron-positron pair creation. The VHE gamma-ray 
absorption, which is energy dependent and increases strongly with redshift, distorts the VHE spectra observed from distant objects. 
The observed energy spectra of the AGNs carry, therefore, an imprint of the EBL. The detection of VHE gamma-ray spectra of 
distant sources (z = 0.11 - 0.54) by current generation Imaging Atmospheric Cherenkov Telescopes (lACTs) enabled to set strong 
upper limits on the EBL density, using certain basic assumptions about blazar physics. In this paper it is studied how the improved 
sensitivity of the Cherenkov Telescope Array (CTA) and its enlarged energy coverage will enlarge our knowledge about the EBL 
and its sources. CTA will deliver a large sample of AGN at different redshifts with detailed measured spectra. In addition, it will 
provide the exciting opportunity to use gamma ray bursts (GRBs) as probes for the EBL density at high redshifts. 
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1. Introduction 

1.1. Cherenkov Telescope Array 

• The proposed Cherenkov Telescope Array (CTA)Q fH lit] is a 
large array of Cherenkov telescopes of different sizes for the de- 
tection of very-high energy gamma-rays. CTA will cover an en- 
ergy range from 10s of GeV up to 100s of TeV with an unprece- 
dented sensitivity - a factor of 10 or more irnprovement over 
current generation instruments such as H.E.S.Sq MAGIcH, and 
VERITASQ- and provide an excellent precision in energy (or- 
der of 10-15%), angular (down to arcmin scale) and temporal 
(down to seconds) resolution fT]. Different layouts of CTA with 
different performance characteristics are currently under inves- 
tigation to balance the financial budget vs. the science return. 
The different CTA layouts and their estimated performance are 
described in detail in a separate paper in this issue [2]. The 
present paper is focused on the projected performance of partic- 
ular CTA layouts in respect to studying the extragalactic back- 
ground light and cosmology. It is demonstrated that CTA has 
a large potential to provide a significant contribution to such 
studies making them one of the central science drivers for the 
instrument. 
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1.2. The extragalactic background light 

During the star and galaxy formation history diffuse extra- 
galactic radiation in the ultraviolet (UV) to far-infrared wave- 
length (FIR) regimes was accumulated. This radiation field, 
commonly referred to as the extragalactic background light 
(EBL), is the second largest background after the Cosmic Mi- 
crowave Background of 2.7 K (CMB) in terms of energy con- 
tained. While the CMB conserves the information on the struc- 
ture of the universe at the moment of decoupling of matter and 
radiation following the Big Bang (at redshift z ^ 1000), the EBL 
is an integral measure of the entire radiative energy released 
dominantly by star formation that has occurred since the de- 
coupling. The EBL energy density from UV to FIR is expected 
to display a two peak structure: a first peak around ~ 1 pim 
from integrated starlight and a second peak at ~ lOOyum from 
infrared emission from dust, which reprocesses the direct emis- 
sion in the UV to optical. 

The precise level of the EBL density is not well known (see 
fs*] for a review). It is difficult to measure directly due to 
strong foregrounds from our solar system and the Galaxy [9J. 
Lower Umits on the EBL density can be derived from inte- 
grated deep galaxy counts at optical and infrared wavelengths 
[e.g. 10, 1 1, 12, 13]. The observation of distant sources of very 
high energy (VHE, typical E > 100 GeV) y-rays using Imaging 
Atmospheric Cherenkov Telescopes (lACT, such as H.E.S.S., 
MAGIC, or VERITAS) provides a unique way to derive indi- 
rect limits on the EBL density. This method will be discussed 
in more detail in the next section. The precision of the EBL con- 
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Figure 1 : Measurements of and limits on the EBL density in comparison to different EBL models/shapes utilized in this study at z = (dashed dotted line: fast 
evolution model from L3J; sohd line: [4J; dashed line: [5J). References for the collection of measurements can be found in [6J and [7J. 



straints set by the lACTs improved remarkably in the last few 



t by the lAL 



years (e.g. 014i l6l ll5ll \ Recently, data in the GeV energy range 
from the Fermi-LAT instrument provided new tools to derive 
constraints from VHE data (e.g. lll6lll7ll ). Contemporaneously 
with the lACT constraints, there has been rapid progress in re- 
solving a significant fraction of this background with the deep 
galaxy counts at infrared wavelengths from the Spitzer satel- 
lite (e.g. ll IM12III8II ) as well as at sub-millimeter wavelengths 
from Herschel or the Submillimeter Common User Bolome- 
ter Array (SCUBA) instrument (e.g. ifioll ). Furthermore, the 
trends of source counts at the faintest magnitudes as measured 
by HST and 8m-class ground-based telescopes (Keck, Subaru, 
VLT) indicate that the EBL contribution from standard galaxies 
is largely resolved [ 10^20] . The current status of direct and in- 
direct EBL measurements are summarized in an accompanying 
review article in this journal and can also be found in |j7|]. 

In total, the collective limits on the EBL between the UV and 
FIR confirm the expected two peak structure, although the ab- 
solute level of the EBL density remains uncertain by a factor 
of two to ten (Fig. [TJ. In addition to this consistent picture, 
several direct measurements in the near IR have also been re- 



ported (e.g. 02111 '). significantly exceeding the expectations from 
source counts (see fS] and f22[] for reviews). If this claimed ex- 
cess of the EBL is real (but see ll23lll4i1 ). it might be attributed 
to emission from more exotic sources like, e.g., the first stars in 
the history of the universe. 

Alternative models do exist, which could, in principle, rec- 
oncile the observations with a high level of the EBL: axion 
like particles (ALPs), which had been introduced to address the 
strong-CP violation problem in particle physics, could lead to 
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Figure 2: EBL attenuation for VHE y-rays for different EBL models and red- 
shifts (z) (FR: 15]; ST: fast evolution model from |3]; GE: |4]). Color code of 
the curves is related to different redshifts (the darker the further away). 



a conversion of VHE gamma-rays to ALPs. This mechanism 
and the related physics for CTA is discussed in a separate paper 
in this special issue fl^. Recently, the secondary gamma-rays 
produced along the line of sight by the interactions of cosmic- 
ray protons with background photons have been invoked to ex- 
plain the VHE spectra of distant sources ll25ll . 
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Figure 3: Energy at which an optical depth of t =1, 2, and 4 is reached vs. 
redshift for the EBL model from |29] (KN, minimal EBL). Also shown are 
approximate energy ranges of CTA (20GeV-100TeV; shaded) and Fermr'-LAT 
(lOOMeV-300 GeV; hatched). 



1.3. VHE y-mys from blazars as a probe of the EBL 

On the way from the source to the observer, VHE y-rays can 
suffer absorption losses by interaction with the low energy pho- 
tons of the EBL via the pair-production mechanism |,26„ ,27ll : 



TvHE + TeBL ■ 



(1) 



VHE y-ray photon fluxes measured on earth (i>obs{Ey) are, 
therefore, altered from the flux emitted at the location of the 
source <I'int(£'y) 



Oobs(£'y,Z) = e 



xOi„,(£,). 



(2) 



The optical depth, T{Ey), encountered by VHE y-rays of en- 
ergy E emitted at redshift z and traveling toward the observer 
can then be calculated by solving the three-fold integral (see 
e.g. ii): 



T(Ey, z) = J d^(z') J dji -!-y^ de n{(E' ,z')cr„(£' ,E'^,ii) 



(3) 



where // = cos 0, n (e) is EBL photon number density and d{(z) 
the distance element (expressed in co-moving coordinates). 

Fig. |2] displays the resulting attenuation due to the EBL for 
VHE y-ray sources at different redshifts and for different EBL 
models (see Fig.[T]). Several features can be identified: 

• At low energies (< 80 GeV) the spectrum is practically 
unattenuated (with the exception of extreme EBL models 
at high redshifts, i.e. the fast evolution model from Stecker 
et al. at z = 2). Since this energy range will be sampled by 
the CTA with high statistical precision, it will be possible 
to measure the unabsorbed spectrum. 

• At energies from 80 GeV to 2 Te V the strength of the atten- 
uation is increasing due to the EBL photons in the optical 
to near-infrared range peak of the EBL density. 



• At energies between 2 to lOTeV a flattening of the atten- 
uation is expected at low redshifts (as seen for z = 0. 1 in 
Fig.|2]for FR and GEN models), due to the ~ /I ' behavior 
of the EBL density in the near to mid-infrared, resulting in 
a constant attenuation. Such a modulation of the EBL at- 
tenuation has been considered as a possible key signature 



for EBL attenuation [e.g. [30|, |3]J, |32|, [16]. Unfortunately 



the intrinsic weakness of the sources combined with the 
sensitivity of the instruments makes it difficult to probe 
such a feature with previous or current generation experi- 
ments (but see llssll '). 



At energies around lOTeV the turnover in the EBL den- 
sity towards the far-infrared peak of the EBL results in a 
strong attenuation, effectively resulting in a cut-off in the 
measured spectra. 



The strength and the position of these features vary with the 
distance of the VHE sources, the assumed EBL model, and the 
overall EBL density. This can also be seen in Fig. [3] which 
shows the energies where the EBL gets optically thick for dif- 
ferent redshifts (i.e. the optical depth reaches unity or higher). 

It should be noted that only distant extragalactic VHE y-ray 
emitters suffer from the absorption by the EBL. For galactic 
sources the attenuation effect from the EBL is negligible up to 
energies of about 100 TeV For higher energies, the absorption 
by the photon field of the CMB starts to be important. Galactic 
diffuse photon fields also produce a weak attenuation signature, 
which might be seen with CTA 13411 . 



The produced e^/e" pairs can initiate electromagnetic cas- 
cades by upscattering low-energy CMB photons to high ener- 
gies; these latter photons would produce, yet again, e^/e~ pairs 
by interacting with EBL photons (see, e.g., llssi [3^). The y- 
rays resulting from the cascade would mostly appear at lower 
energies than the ones accessible by CTA but they would mod- 
ify the overall energy spectra of the sources, which is impor- 
tant for their modeling. If the strength of intergalactic magnetic 
fields and interaction length are low enough, the cascade com- 
ponent would appear to arrive from the direction of the source. 
Higher magnetic fields would isotropize the cascade compo- 
nent. 

Since the EBL attenuation alters the VHE y-ray spectra of all 
distant sources, understanding the EBL's properties are crucial 
for studying the intrinsic source physics. On the other hand, 
VHE y-rays from distant sources can be used as probes for 
the intervening photon fields. This requires some assumptions 
about the intrinsic spectra, which, e.g., can be derived from 
nearby sources or unattenuated parts of the spectrum. In ad- 
dition, EBL attenuation should be a global phenomenon affect- 
ing the spectra of all sources in a consistent way. In particular, 
sources at the same distance should show the same attenuation 
imprint from the EBL. 
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2. The CTA EBL & cosmology physics case 

2.1. Disentangling intrinsic blazar spectra from EBL effects: 
precision measurement of the EBL density 

There is a fundamental difficulty to distinguish between in- 
trinsic blazar effects (e.g. intrinsic cut-offs or possible multi- 
zone emission regions) and the absorption effect caused by the 
EBL. This ambiguity is also the main reason why, with al- 
ready more than 40 extragalactic sources measured by H.E.S.S., 
MAGIC and VERITAS, only upper limits on EBL density have 
been derived yet (as in, e.g., 1 14]). The number of sources de- 
tected at VHE will increase further in the next 5 years, in partic- 
ular owing to the upgrades of the existing facilities aiming at en- 
hanced sensitivity and lower the energy thresholds: MAGIC II, 
H.E.S.S. II and VERITAS II. However, the real boost in the field 
will be achieved with CTA only, when it will become possi- 
ble to measure simultaneously and with high precision both the 
unattenuated (usually at energies below 100 GeV) and attenu- 
ated (at higher energies) parts of the AGN spectra. Such mea- 
surements will allow to disentangle between intrinsic features, 
which can vary from source to source and even for the same 
source depending on its flux state, and the EBL induced signa- 
ture, which solely depends on the redshift of the source. Also, 
precise measurements at energies > 10-20 GeV (the foreseen 
energy threshold of CTA) in combination of multiwavelength 
data would ease and improve the modeling accuracy of AGNs 

Disentangling the intrinsic AGN features and the EBL in- 
duced attenuation signature in the measured AGN spectra will 
unambiguously lead to a precise measurement of the EBL den- 
sity at a redshift of z = 0. The precision of the measure- 
ment will depend on the richness of the AGN sample avail- 
able, in particular on the number of sources at different red- 
shifts to enable a robust modeling of the intrinsic features of 
the sources. The power of CTA to simultaneously measure the 
intrinsic and the EBL attenuated part of an energy spectrum is 
investigated in Section [3.2.1l utilizing the measured VHE of the 
blazar PKS 2155-304 as a template. 

In individual sources an intrinsic break or cut-off cannot be 
excluded. However, with CTA it is estimated (see |38]) that 
over a period of 5 years conservatively 100 new AGN sources 
could be expected to be discovered. This gives a sizable source 
sample to check for global features/trends in the spectrum as 
a function of the expected opacity (see, e.g., IitIi for a study 
on the current VHE AGN sample), or redshift (see, e.g., 1139 1 
for a simulated AGN sample) allowing to exclude biases due 
to source specific features. Combined with statistical studies of 
a large number of objects, precision measurements point to a 
unique and powerful method to probe the EBL that is difficult 
to achieve with present-day instruments. 

Owing to the dependence of the pair-production cross section 
on the energy of the participating photons, distinct AGN types 
are crucially important in the indirect EBL measurements at dif- 
ferent wavelengths. For measuring the EBL in UV, the sources 
must be distant to experience a measurable absorption in UV 
and they must be sufficiently bright to provide enough flux for a 
CTA measurement. Moreover, their intrinsic spectra must reach 



energies up to 300-500 GeV in order to connect from UV to the 
optical regime. Therefore, distant bright AGNs, mostly FSRQs 
(e.g. 3C279 |15] and PKS 1222 |40]), and maybe GRBs if ex- 
tending to high enough energies are the sources with the best 
potential to probe the EBL in UV. For measuring the EBL in 
optical to mid IR regimes, AGNs with hard spectra extending 
to few TeV are the ones providing the best results. The red- 
shift of the sources must not be too large (z<0.5) in order to 
avoid a total absorption of the flux at energies above 1 TeV. 
This is a guaranteed science case, since good candidates with 
hard energy spectra already exist: lES 0229-1-200 RGB 
J0710+591 lEl, lES 1101-232 (a^. A similar case can be 
made for the far IR regime. There, few close by sources (e.g. 
Mrk 421 |44], Mrk 501 H, M87 iS], IC310 il) are the 
best ones provided they produce intrinsic spectra up to energies 
of at least E>20TeV. 

The promising technique of the indirect EBL measurements 
can be combined with measurements from deep source counts 
at optical to IR wavelengths. Through this, the contribution 
to the EBL from faint unresolved sources or any exotic pop- 
ulations can be disentangled from relatively bright and known 
populations. 

2.2. Star formation history and the early universe 

The cosmic star formation rate density (SFRD), one of the 
fundamental quantities of cosmology, strongly evolves over the 
redshift range z = to z = 2 with a peak expected at redshifts 
between 1 and 2. While the SFRD up to z < 1 is reasonable well 
measured from different tracers, at higher redshifts the uncer- 
tainties are large. Since stellar and dust emissions are expected 
to be the dominant contributors to the EBL, a precise measure- 
ment of the EBL density can provide strong constraints on the 
SFRD iRsll . The science output is guaranteed since the SFRD 
is a highly debated topic of the modern cosmology. It is esti- 
mated that the precision of the SFRD determination, which can 
be achieved using VHE spectra, will be on a level of 20-30% 
(based on iH]). 

The end of the dark ages of the universe - the epoch of reion- 
ization - is a topic of great interest fsO*, 'sT]. Not (yet) accessi- 
ble via direct observations, most of our knowledge comes from 
simulations and from integral observables like, e.g., the cosmic 
microwave background. The first (Population III) and second 
(Population II) generations of stars are natural candidates for 
the sources of reionization. If the first stars are hot and mas- 



sive, as predicted by simulations 11521 15311 . their UV photons 
emitted at z > 5 would be redshifted to the near-infrared regime 
and could leave a unique signature on the EBL density 15411 . If 
the EBL contribution from lower redshift galaxies is sufficiently 
well known (e.g., as derived from source counts) upper limits 
on the EBL density can be used to probe the properties of early 
stars and galaxies [55']. Combining detailed model calculations 
with redshift dependent EBL density measurements could en- 
able to probe the reionization history of our universe. Detec- 
tion of spectral cutoffs in sources at very high redshifts (z > 5) 
would directly probe the evolving UV background at these red- 
shifts, providing invaluable insight into the cosmic reionization 
epoch. 
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While blazars may be difficult to detect much beyond z « 2, 
this may be feasible for GRBs. GRBs are the most luminous 
and distant gamma-ray sources known in the Universe, typi- 
cally arising at redshifts z ~ 1-4, which correspond to the 
peak epoch of cosmic star formation activity. Furthermore, they 
are known to occur at least up to z ~ 8 - 9 1561 571. well into 



the cosmic reionization era, and possibly even beyond, out to 
the very first epochs of star formation in the Universe. The 
recent detections by FermZ/LAT of dozens of GRBs including 
GRB 080916C at z = 4.35 158J clearly demonstrate that at least 
some GRBs have luminous emission extending to few tens of 
GeV lf59ll60l] . Moreover, the duration of the multi-GeV emis- 
sion is seen to last up to several thousand seconds, in agreement 
with reasonable theoretica l exp ectations for the high-energy af- 
terglow emission see e.g. [60]. Thanks to CTA's much greater 
effective area compared to Ferm;/LAT (by about a factor of 10^ 
at 30 GeV if using CTA array layout B) together with its rapid 
slewing capability (180 degrees within 20 seconds), there are 
good prospects for high precision measurements of the GeV- 
TeV spectra of GRBs at z > 1, as well as the detection of GRBs 
at z ^ 4 and beyond, a regime that cannot be explored with 
AGNs. 

2.3. Particle physics and the unknown 

Data collected so far by H.E.S.S., MAGIC and VERITAS on 
absorption for AGNs have constrained the maximum EBL den- 
sity to a level close to the lower limits from integrated galaxy 
counts. Once CTA goes online and significantly improves the 
indirect EBL measurements, several scenarios are possible. In 
the first scenario the joint efforts of direct and indirect measure- 
ments converge and the EBL density is resolved. If the EBL 
density is well measured, EBL absorption in VHE spectra can 
be used as a distance indicator, somewhat analogues to what is 
done for Supernovae o f typ e 1 A but with very different system- 
atic uncertainties (see f49']). The evolution of such an attenua- 
tion signature with redshift would allow us to put constraints on 
cosmological models and measure the Hubble parameter and 
cosmological densities. The proposed method would provide 
independent and complementary constraints on the cosmolog- 
ical parameters with different systematic uncertainties to what 
projects as LSST |61] or PLANCK I'iS] are planning to reach. 
An independent distance measurement of AGNs using the EBL 
absorption feature in their VHE spectra is also important for 
getting redshifts for many BL Lacs, since only 50% of all BL 
Lacs have reliable redshift measurements. 

In a second scenario, the EBL level determined from the 
VHE observations will be remain to be a factor of two to three 
higher than the level resolved through the galaxy counts. This 
would point to a large population of yet hidden sources, maybe 
also to a truly diffuse background such as due to decaying big- 
bang relict particles. In case the discrepancy between the lower 
and upper limits will be restricted to a narrow wavelength range, 
the search for the missing background photons will be easier. 

In an alternative scenario, the upper limits on EBL from in- 
direct measurements will contradict lower limits from galaxy 
counts. This opens up a completely new window in particle 
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Figure 5: Example of simulated spectra for dilferent EBL densities. The 
base spectrum assumed is the quiescence state spectrum of PKS 2155-304 
(z = 0.116), EBL model is from |5], used with different scalings. Shown are: 
the measured spectrum (grey markers), the simulated spectra for different level 
of the EBL density (black markers) and the corresponding assumed intrinsic 
spectra (black lines), the source spectrum in the GeV energy range as measured 
by Fermj-LAT |65] (not simultaneous, purple butterfly), and the energy ranges, 
which are used to determine the slope of the simulated spectrum at low (light 
red) and high (blue) energies. For details on the data and methods see i66ll . 



and astroparticle physics since one of the two following op- 
tions must be true: either the intrinsic source physics is funda- 
mentally different than the models predict or the properties of 
gamma ray propagation through space are fundamentally differ- 
ent than we know today, which can become a major challenge 
for the modem physics (see e.g. (63l 24ll ). 



3. Investigating the CTA response 

3.1. Simulations 

For the studies presented here the CTA physics response sim- 
ulation tools presented in [64] are used. In short, a source in- 
trinsic spectrum is assumed in each case and then folded with a 
given CTA response. The source is located at the center of the 
telescopes' field of view and the background (especially impor- 
tant at low energies) is assumed to be measured with a 5-times 
better statistics than the exposure of the source region. Only 
significant, i.e. larger than 3 standard deviations, reconstructed 
flux points have been considered in the analysis. Furthermore, 
the systematic uncertainty is considered to be 1 % of the back- 
ground level, and we require the signal to be at least 3 times 
higher than that. The performance files for 20deg zenith angle 
(fixed) are used and no correction for source motion along the 
sky is made. 

An example of such a simulated spectrum is shown in Fig.|3] 
where the measured spectrum of the blazar PKS 2155-304 is 
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Figure 4: Results from the power law fit to the simulated spectrum in the low and high energy band using the PKS 2155-304 quiescence spectrum, the EBL model 
from |5], and CTA aiTay layout B (20h exposure). Left: Attenuation of the VHE y-ray flux resulting from the scaled EBL models (lines). The colored boxes in the 
background mark the fit range for the power law at low (red) and high (blue) energies. Right: Spectral index T from the fit of a power law to the low (red) and high 
(blue) energy range versus the scaling factor of the EBL model. The color shaded bands denote the error on the spectral index from the fit (RMS of the spectral 
index distribution). Black crosses mark the intrinsic spectral index that has been utihzed. For discussion see main text. 



used together with different assumptions about the EBL density 
(here a scaled version of the EBL model by [5]) to calculate 
the spectrum as it would have been detected by CTA in a 20 h 
exposure (CTA array layout: B) (more details can be found in 

m). 

3.2. Results 

3.2.1. Disentangling intrinsic blazar spectra and EBL effects. 

One of the main features of CTA will be a high sensitivity 
in the energy range between 20 and lOOGeV, an energy range 
which holds the possibility to directly sample parts of the en- 
ergy spectrum of a source, which are not affected by the EBL 
attenuation. With the difference between the measured spec- 
trum in the unabsorbed part of the VHE spectrum and in the 
absorbed part of the spectrum, especially if studied for several 
sources with good statistics, the strength of the near-IR EBL 
can be derived. Here, exemplarily, the case of PKS 2155-304 is 
studied following the approach described in ll66ll . 

For different levels of the EBL density, generated by scal- 
ing the EBL model from |5], the intrinsic source spectrum is 
reconstructed and used to simulate the corresponding spectrum 
as measured by CTA (Fig.|5]l. The resulting attenuations for the 
scaled EBL models are shown in Fig.|4]left panel. In the CTA 
simulated spectrum two distinct regions in energy can be iden- 
tified: (i) the unabsorbed region, roughly from 40 to lOOGeV 
for the redshiftofPKS 2155-304 of z = 0.116 (red band), where 
no or little EBL attenuation is present and the intrinsic spectrum 
can be observed directly, and the (ii) absorbed region at energies 
> lOOGeV, where the EBL attenuation can be measured (blue 
band). In each of the two bands the spectrum is characterized 
by a power law dN/dE ~ E^^ and the derived photon indices 



r, shown in Fig.|4]right panel, are then investigated. The fit in 
the low energy band reproduces very well the assumed intrinsic 
spectral indices (red dots vs black crosses in Fig. |4]right panel). 

Fig. |6]displays the strength of the spectral break between the 
two power laws Flow - Fhich for different CTA array layouts, 
with the error bands calculated via error propagation. For array 
layout B (best low energy performance) an EBL density a step 
in scale factor of 0.1 corresponds to one to two standard de- 
viations difference in the strength of the break. Array layout I 
provides a comparable but slightly worse performance while the 
decreased low energy sensitivity of array layout C makes stud- 
ies of this type difficult (resolution factor two to three worse and 
the intrinsic index cannot be reconstructed). Note that a scaling 
of 0.1 corresponds to an EBL density of ~1 nWm^^sr ' at 2/zm 
which is of the same order as the error on the lower limits from 
integrated source counts at this wavelength. Since EBL attenu- 
ation is a global, redshift dependent phenomenon, in general, a 
combined fit to all available VHE data should be used to derive 
limits on the EBL density [see e.g.|^. 

While for many of the up to now detected extragalactic VHE 
sources the energy spectrum in the VHE range is well described 
by the extrapolation of the spectrum in the HE range folded 
with an EBL attenuation some of the sources show an intrin- 
sic break between the two energy regimes (e.g. PKS 2155-304) 
1 6711 . For individual sources, such an intrinsic break could lead 
to a wrongly determined EBL density when blindly applying 
the methods outlined above. The solution, again, lies in uti- 
lizing a sample of source: while intrinsic features are different 
between the sources the EBL attenuation signature can be pre- 
cisely determined for each of them enabling one to disentangle 
the two. 
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EBL density scale factoi EBL density scale factor EBL density scale factor 

Figure 6: Spectral break between the power law fits in the low and high energy bands for different CTA array layouts (spectrum: PKS 2155-304 quiescence; EBL 
model: 0]; exposure: 20h). Circles mark the reconstructed spectral break derived from the simulated data while crosses mark the spectral break between the 
simulated spectral index and the reconstructed index in the high energy band. The error bands are deiived from error propagation. For more details on the anlysis 
see L66i1 . 



3.2.2. Distant sources, GRBs 

Due to the low energy threshold of CTA and its high sensitiv- 
ity CTA will have the opportunity to detect even distant GRBs. 
Whereas a discussion on GRB physics is beyond the scope of 
this article, GRBs start playing a crucial role in understanding 
EBL density at high redshifts [68]. A common expectation is 
that only GRBs will provide high enough y-ray luminosity to 
detect a source located at high redshifts (z > 2). Also AGN in 
extreme flaring states could be detected at z > 2 (see Figure 6 in 
AGN paper, same issue). As mentioned before, y-ray sources 
beyond a redshift of two can be used as beacons probing the 
star formation rate density and, in particular, its peak position. 

The power of CTA to detect distant GRBs can be demon- 
strated by utilizing the ferm/-LAT detection of the extraordi- 
nary GRB 080916C, which took place at a redshift of z = 4.3 
lISSll . In Fig.|2]we compare simulated spectra as they were mea- 
sured by CTA of the GRB. In the left panel, we assume the 
instrinsic spectrum as it was measured by Fenni-LAI for the 
time-span To = 55 - 100 s (45 s, interval "e" as defined in 
ifssll l The result is that a clear detection can be made, here 
on an example of CTA array layout B (see [2] for details on 
individual array configurations. For this array layout, it would 
also have been possible to derive a detailed spectrum. Com- 
paring performance of different array layouts for this physics 
case, we see only marginal improvement using the array lay- 
out B (low threshold) instead of £ or / (balanced arrays): the 
resulting error bars are only few percent better and the recon- 
structed energy range is identical. However, for the array layout 
C (high threshold), no detection of this GRB can be achieved 
and, consequently, no spectrum can be derived. Since GRBs 
are short lived events, a high repositioning speed of CTA is 
required. In the right panel of Fig. |7]we see a simulation of 
the GRB080916C spectrum for the interval "d" L58jfl The in- 



The assumed 



source flux is dN/dE = 1.4 X 
'TeV^', whereas the index -1.85 is derived from 



10"''(£/TeV)-'*'^cr 
a fit to the LAT-only data (Figure S4 in the Supplement to 
communication from the Fermi-LAT team). 



and private 



trinsic spectrum is harder and the flux level is higher than for 
the interval "e", which results in clear detection even for an in- 
tegration time of 20s. Results for four different EBL models 
Kneiske et al. "best fit" |690 (red), Finke et al. S (green), 
Gilmore et al. |71] (magenta), and Inoue et al. (private commu- 
nication, blue) are presented, which show clear difference and 
permits distinguishing between the models. We note that catch- 
ing a GRB at a time To < 55 s after the onset of the GRB, which 
would correspond to the interval "d" in case of GRB080916C, 
is challenging even if CTA would be able to reposition to any 
pace in the sky within 20 s. It is, however, withing the reach 
of the instrument. Observations at Tq > 55, corresponding to 
the interval "e" in case of GRB080916C, are easier to achieve 
and can be considered as a secured case. It is important to men- 
tion that GRB080916C was the brightest GRB observed by the 
Fermi/LAT so far. The probability that GRBs with sufficiently 
high ffux will be observed by CTA within its life time has been 
addressed by several authors 

M III finding that CTA, if hav- 
ing a low energy threshold of 20-30GeV would be able to de- 
tect 0.1-0.2 GRBs per year during the prompt phase and about 
1 per year in the afterglow phase. 

As discussed by a number of authors ll59[ l60[ [tsIi . UV radi- 
ation with sufficient intensities to cause the reionization of the 
intergalactic medium (IGM) are also likely to induce apprecia- 
ble gamma-ray absorption in sources at z > 6 at observed ener- 
gies in the multi-GeV range, with a potentially important con- 
tribution from Pop III stars. Measurements of these effects can 
thus provide important cross-checks of current models of cos- 
mic reionization, a unique and invaluable probe of the evolving 
UV EBL during the era of early star formation, as well as a test 
for the existence of the yet hypothetical Pop III stars. For more 
details we refer to the GRB paper in this issue. 

3.2.3. Precision measurement of today's EBL 

The measured energy spectra of AGNs in the energy range 
between lOOGeV and few TeV follow usually a smooth shape. 
For most of the measured sources, a simple power law fit is suf- 
ficient to describe the available data well, whereas for sources 




Figure 7: Simulated energy spectrum of GRB 080916C (z = 4.3) if measured witli CTA. Left panel: Using the CTA array layout B (low energy threshold) and 
assummg EBL from (vl. The intrinsic spectrum is assumed to follow 2.9 X 10"' X (£/TeV)"2.i6|--pgY-i 

cm s 1 and the duration of the measurement is 45 s 
(To=55-100s, interval "e") as measured with Fermi/LPCT |58]. A clear detection can be made and spectral shape can be measured from 30 to 100 GeV. The result 
is basically identical if using CTA array layouts E or /. For the array layout C (high energy threshold), no detection of this GRB can be made. Right panel: Using 
the CTA array layout E and exposure time 20 sec for the interval "d" ISSJ, i.e. the assumed source MxiidNldE = 1.4 X 10"'(£/TeV)-''*5cm-2s-lTeV-'. Here 
we compare the EBL models of Kneiske et al. "best fif [69] (red), Finke et al. iTOh (green), Gilmore et al. ItiIi (magenta), and Inoue et al. (private communication, 
blue). For the array layout C, again, no detection of this GRB can be made. 



in a flare state (like the flai-e of PKS 2155-304 in 2006) or 
with a generally high emission state (like Mkn421), either a 
curved power law or a power law with a cut-off are successfully 
used. The curved power law (also known in the literature as the 
double-log parabola) is expected to describe the spectra well at 
energies close to the position of the Inverse-Comptonpeak. The 
power law with a cut-off instead is the expected behavior of a 
source which does not provide necessary conditions for acceler- 
ation of charged particles to sufficiently high energies. All sce- 
narios, assuming the one-zone model, which seems adequate 
to describe most VHE spectra of AGNs rather well, do have 
one common feature: the measured spectra can be described by 
smooth functions, i.e., no features, wiggles or pile-ups are ex- 
pected, especially after de-convolving the spectra for the effect 
of the EBL absorption. This property can, therefore, be used to 
distinguish between different overall EBL levels in the optical 
to infrared regime: whereas the "correct" EBL model and level 
will produce a smooth intrinsic AGN spectrum, an "incorrect" 
EBL level would result in a signature (in form of well defined 
wiggles) in the reconstructed intrinsic spectrum SS. Pos- 
sible spectral features in the intrinsic spectra could complicate 
the analysis, though other methods like investigating the spec- 
tral properties of variable sources on short time-scales or us- 
ing multiwavelength based modeling might help to disentangle 
source intrinsic and EBL effects (see Sec 2.1). It is also impor- 
tant to observe the same effect on spectra of different sources to 
become independent of possible intrinsic source features. An- 
other caveat is that the form of the wiggles will have some 
dependency on the exact spectral shape of the "correct" EBL 
model, which is not known a-priori. Though the state-of-the- 
art EBL models agree within 10% on the expected EBL shape 
(with the only exception being the model of |3]), a possibility of 
a significantly different shape, e.g., due to a contribution from 
early stars, cannot be excluded. 



The strength of the method is that the EBL signatures in the 
reconstructed AGN spectra will not only be visible (measur- 
able) in the case where the assumed EBL level is higher than 
the real one, but also in case the assumed EBL level is lower 
than the real one. It is, therefore, the first indirect method to 
really measure the EBL density at z-0. 

Simulation example. The method is illustrated in Fig. |8]for the 
VHE spectrum of Mkn501. The assumed spectral shape and 
the flux level of the intrinsic spectrum are adapted to the average 
flux measured by HEGRA |76] during the outburst of the source 
in 1997 (llOhobservation time): the original data are shown by 
grey open squares in the upper panel of tiie figure. The simu- 
lated CTA spectrum calculated using the 'correct EBL' (model 
from [^) is shown in red, whereas the assumed intrinsic spec- 
trum of Mkn 501 is shown by the solid grey line and the recon- 
structed intrinsic spectrum is shown by the blue filled circles. 
The effect of the mis-reconstruction of the intrinsic spectrum is 
shown for the example of an EBL scaled by a factor of 1.3: the 
reconstructed intrinsic spectrum (green filled triangles) clearly 
shows wiggles in the fit range. The effect of the wiggles is more 
visible in the lower panel of the figure where the residuals to the 
best fit function are shown. The wiggles are quantified by a fit 
in the energy range between 100 GeV and 7 TeV, well before a 
possible pile-up in the spectrum arises, which can be addition- 
ally used to rule out an EBL realization lE8[ l6ll. The choice of 
the fit range is made in order not to bias the result by the level of 
the EBL above 10 /vm, to which the VHE spectra are very sen- 
sitive due to a super exponential dependency of the attenuation 
with the wavelength in that range. Using the correct EBL level 
to reconstruct the intrinsic spectrum, the intrinsic spectrum is 
well described by a smooth function (Fig.|8] middle panel). In- 
stead, when using a "wrong" scaled EBL density characteristic 
deviations (wiggles) from a smooth function are visible (Fig.[8l 
lower panel). The wiggles are quantified by a reduced j^'^ value 
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Figure 8: Search for signatures at mid-VHE in the Mkn 501 spectrum. Top panel: The MJmSOl spectral energy distribution (SED) at very high energies. Shown 
are: measured spectrum by HEGRA (grey open squares), simulated CTA spectrum assuming the array layout / (red filled squares), assumed intrinsic spectrum of 
Mkn 501 (grey solid line), reconstructed intrinsic Mkn 501 spectrum (blue filled circles) and reconstructed intrinsic spectrum of Mkn 501 assuming an EBL scaling 
factor of 1.3 (green triangles). The Crab Nebula spectrum (grey dashed fine) is shown for compaiison. Middle panel: Residuals between the best fit to the SED 
in case of CTA measurement and the spectral points from the intrinsic spectrum using the correct EBL density. Bottom panel: Residuals between the best fit to 
the SED for a CTA measurement and the spectral points from the intrinsic spectrum reconstructed using the scaled EBL model. A clear and significant signature 
(wiggles) in the residuals is visible, which is quantified by a low probability of the fit. 



of 61.4/12, corresponding to a fit probability of 1 .3 • 10 The 
small fit probability for the assumed scaled EBL level implies 
(under used assumptions) significant presence of the unphysi- 
cal wiggles and, therefore, the exclusion of that particular EBL 
realization. 

Results of the analysis. This procedure is repeated 1000 times 
for every scaled EBL density in order to achieve a solid statis- 
tical mean and a 1 cr (68%) coverage of the reduced values. 
The results are shown in Fig. |9]for two different CTA exposures 
(20 and 50 hours) plotting the mean reduced values includ- 
ing 68% error bars versus the EBL scaling factor As a result 
one can see a parabola-like curve with the minimum at EBL 
scaling factor of 1.0, which is via construction the correct EBL 
model. The expected mean reduced values are shown by the 
blue filled squares and the green filled circles for the exposures 
of 20 hr and 50 hr, respectively. The corresponding fit probabil- 
ities of P=50%, 1% and 0.01% are shown by the dotted, short- 
dashed and long-dashed lines, respectively. An EBL scaling is 
considered to be excluded when the mean reduced value in- 
cluding its 68% error exceeds the x^ value for P=0.01%. For 
the case constructed here this means that in case of 50 hours 
observation, the EBL scalings below 0.70 and above 1.30 are 



excluded. In the case of 20 hours observation, the EBL scalings 
factors above 1.40 and below 0.6 are excluded. 

When comparing performance of different CTA array layouts 
for this particular study, no significant differences between the 
results was found for the array layouts B, I, and C. The rea- 
son for low discrepancies is that the study mainly relies on the 
energy resolution of the array in the range 200 GeV up to few 
TeV and per construction this is the regime where all CTA ar- 
rays have similar performances. 

4. Conclusions, recommendations, & caveats 

CTA will be a superb machine for probing the EBL. Key fea- 
tures are the high sensitivity, providing a large and rich sample 
of sources for EBL studies, and the huge energy coverage, en- 
abling to simultaneously measure the EBL attenuation and the 
un-attenuated intrinsic spectrum. CTA will deliver definite an- 
swers to some of the most pressing topics in this field, including 
a precise measurement of the EBL density and it evolution up to 
redshifts of one and beyond, providing new insights into struc- 
ture formation and galaxy evolution. In addition, it will open up 
the field of EBL studies towards more fundamental physics, al- 
lowing to probe for effects of fundamental particles like hidden 
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• For low energy studies of rapid transient phenomena like 
GRBs, no significant difference between the performance 
of arrays B and / could be found. Due to the short expo- 
sure, the background rejection does not play an important 
role, equalizing the performance of the two types of arrays 
(see Sec. l3.2.2] l. Instead, an array of type (3) (e.g. C) with 
a significantly higher energy threshold does not allow for a 
detection of the phenomena and is, therefore, unsuited for 
such studies. 

• For the precision measurements of the z = EBL density 
as discussed in Sec. 13.2.31 all tested arrays show very sim- 
ilar performance. The EBL can be resolved down to 20- 
30% thanks to the superb spectral resolution of the planned 
CTA. However, in order to be sensitive to the differences 
between the state-of-the-art EBL models, an energy reso- 
lution of 5-10% in the energy range between 100 GeV and 
few TeV is required. 



Figure 9: Quantitative results from the search for EBL signatures in the mid- 
VHE using the energy spectrum of MknSOl. The results for the array layout 
/ are shown: the reduced values from the fits to the reconstructed spectra 
of MknSOl as a function of EBL scaling factor Blue filled squares and green 
filled circles show the expected result for 20 and 50 hours of CTA observa- 
tions, respectively. The black horizontal lines correspond to fit probabilities as 
labeled. 

photons or axions (see [^24"]), and give constraints on cosmol- 
ogy- 

The currently proposed layouts come in three different fla- 
vors: (1) emphasis on low energy sensitivity and lower thresh- 
old with worse high energy performance (array layout B), and 
(2) increased high energy performance with worse acceptance 
at low energies (array layout C), and (3) overall balanced sensi- 
tivity with emphasis on the core energy range between 100 GeV 
and 10 TeV (e.g. array layout E or /). Each of the setups have 
their advantages for certain topics of the EBL studies: arrays of 
type (1) are thought to be better suited for studies of UV radi- 
ation and early universe using GRBs and distant AGNs where 
the most observational emphasis is put on energies below 100 
GeV; array layouts of type (2) provide the opportunity to study 
in great detail the mid and far IR with sources like M 87 and 
lES 0229-1-200 where VHE gamma ray emission extends in en- 
ergy up to 10s of TeV; and while array layouts of type (3) seem 
to deliver the best compromise of the previous two types. 

Several CTA array layouts have been investigated in detail 
for their performance for different EBL and cosmology physics 
cases. The results are: 

• For low energy studies of steady objects (i.e. implying 
long exposure), array of type B is superior in its perfor- 
mance thanks to its large collection area and good back- 
ground rejection at energies below 100 GeV. Due to their 
decreased sensitivity at lower energies, the capability of 
sampling the un-attenuated part of the spectrum with array 
layouts / and E is slightly worse (see Sec l3.2.l] i. Array 
layout C shows significant disadvantages for these type of 
studies. 
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